Glassy Vortex State in a Two-Dimensional Disordered XY-Model 
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The two-dimensional XY-model with random phase-shifts on bonds is studied. The analysis is 
based on a renormahzation group for the replicated system. The model is shown to have an ordered 
phase with quasi long-range order. This ordered phase consists of a glass-like region at lower 
temperatures and of a non-glassy region at higher temperatures. The transition from the disordered 
phase into the ordered phase is not reentrant and is of a new universality class at zero temperature. 
In contrast to previous approaches the disorder strength is found to be renormalized to larger values. 
Several correlation functions are calculated for the ordered phase. They allow to identify not only 
the transition into the glassy phase but also an additional crossover line, where the disconnected 
vortex correlation changes its behavior on large scales non-analytically. The renormahzation group 
approach yields the glassy features without a breaking of replica symmetry. 
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I. INTRODUCTION 



The two-dimensional XY-model with random phase- 
shift, defined in Eq. (|l|) below, captures the physics 
of a variety of physical systems. Among them are 
XY-magnets with non-magnetic impurities, which are 
coupled to pthe XY-spins via the Dzyaloshinskii-Moriya 
interaction,!!! Josephson-junction arrays with geometric 
disorder,H and crystals with quenched impurities.tl Dis- 
ordered XY-models are related to even more physical 
systems, like impure superconductors. It is in partic- 
ular the interest in vortex-glass phases in type-II su- 
perconductors, which motivates an analysis of possi- 
ble glassy features in the paradigmatic two-dimensional 
bond-disordered system. 

The most fundamental question is, whether this model 
has an ordered phase. In the absence of disorder, there 
exists at low temperatures a phase with quasi long-range 
order. The famous Kosterlitz-Thouless transition sep- 
arates it frora-pthe high-temperature phase with short 
ranged order.Ql3 The analysis by Cardy and Ostlund 
and more explicitly by Rubinstein, Shraiman and NelsonE 
predicted that the phase transition should be reentrar± 
in the presence of disorder. Later Korshunov suggested,!!! 
that the ordered phase might ba completely destroyed by 
disorder. Ozeki and NishimoriE! showed for models with 
gauge-invariant disorder distributions, that the phase 
transition cannot be reentrant, leaving open|Whether the 
ordered phase does exist or not. ExperimentsO on Joseph- 
son junction arrays and simulationalj were in favor of the 
existence of an ordered phase without reentrance. 

In this controversial situation Nattermann,_Scheidl, 
Korshunov, and Lillil as well as Cha and Fertiglj recon- 
sidered only recently the problem and found that the 
ordered phase does exist and its boundary does not have 
reentrant shape. The earlier observation of reentrance 
was attributed to an overestimation of vortex fluctua- 
tions at low temperatures. 

The present study extends Refs. |l^ and ^ which fo- 



cused on the absence of reentrance. The purpose of this 
article is to examine in more detail the ordered phase, 
which is shown to be composed of different regions. They 
are distinguished by the behavior of vortex correlation 
functions on large length scales, which are calculated 
quantitatively. A low-temperature regime is found to ex- 
hibit glassy behavior in its correlations. 

In addition to the previous studies we find that the 
strength of disorder is increased under renormahzation. 
This effect does not destroy the ordered phase but it is 
crucial for the universality type of the transition at low 
temperatures. 

In the following a generalized self-consistent screening 
approach is applied to the replicated system—Tie main 
improvement compared to earlier approachesEj!ji3 on the 
same basis is achieved by including all relevant fiuctu- 
ations and by taking special care when the number of 
replicas is sent to zero. 

We use the replica formalism mainly for two reasons: i) 
The calculation of correlations becomes very convenient: 
one may use an expansion techuictue (small fugacity and 
density of vortices) which failstil^llj in the unreplicated 
system, ii) We show that in the present model there is no 
need for replica symmetry breaking, even at lowest tem- 
peratures. In particular, the reentrance disappears and 
the glassiness occurs without a breaking of replica sym- 
metry. This might astonish, since for various disordered 
systems reentrance has been explained as an artifact of a 
replica-symmetric (exact mean-field or approximate vari- 
ational) approach, which has been overcome by breaking 
the replica symmetry.tl! 

The paper is organized as follows: In Sec. II we set 
up our model and map it onto an effective Coulomb gas. 
We recapitulate the breakdown of fugacity expansion in 
the presence of disorder. In Sec. Ill we set up the renor- 
mahzation group within a self-consistent formulation of 
screening. The consistency of this approach in the orig- 
inal and replicated system is demonstrated. Sec. IV is 
devoted to the evaluation of the self-consistency and the 
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proper performance of the rephca hmit. The results of 
the renormalization group treatment are worked out in 
Sec. V, which are discussed and compared to the results 
of previous work in Sec. VI. 



II. THE MODEL 

The XY-model with random phase-shifts is given by 
the reduced Hamiltonian 

n^-K COs(6lr -6lr' -^<r,r'>)- (1) 

<r,r'> 

It refers to a square lattice with unit spacing in two di- 
mensions. An XY-spin with angle Or S ["""j ""i is at- 
tached to every lattice site (position r). The reduced 
spin coupling between nearest neighbors < r, r' > reads 
K = J/T. The interaction involves quenched random 
phase-shifts with variance ^<i.r'> ~ which are un- 
correlated on different bonds. Thermal fluctuations are 
weighted according to the partition sum 




This XY-model can be mapped approximately onto an 
effective Coulomb gas. SipeeJjhis procedure is well known 
in the absence of disordeili3ll3 and can be performed sim- 
ilarly in the presence of disorder, we recall only briefly 
the main manipulations leading to the effective model. 
In order to simplify the analysis of model ^ , one can re- 
place approsipiately the cosine- interaction by the Villain- 
interaction.ll3 Thereby an additional variable, the vortex 
density N is introduced. The decisive advantage of this 
approximation is, that the Hamiltonian becomes bilinear 
in the angles. All configurations of the original angles can 
be expressed in terms of spin-waves and vortices. Just as 
in the pure case, spin-waves and vortices decouple en- 
ergetically. Since spin-waves are simply harmonic, all 
non-trivial physics arises from vortices with an effective 
Hamiltonian 

= y |A^k - Qkl' (3a) 

= S^(iVR-0R)2- (3b) 
R 

- 7ri^ln|R-R'|(iVR-QR)(iVR, -QrO, 

R/R' 

in Fourier and real space, and the partition sum 

{N} 

Here TVr is the integer-valued vortex density, which is 
subject to thermal fluctuations. On the square lattice TVr 



is associated with plaquettes (dual sites), which are la- 
beled by R in contrast to sites r. Thermal vortices N cou- 
ple to a background of quenched vortices Qr = ^ V x A. 
The last expression means, that the plaquette variable 
Qr is given by the rotation of variable A on the sur- 
rounding bonds. From this relation one immediately de- 
rives, that the quenched vortices are Gaussian random 
variables with variance Q-kQk = if^^^- ^ single phase- 
shift A<r.r'> creates a dipole of quenched vortices situ- 
ated in the plaquettes which have the bond < r, r' > as 
side. Therefore quenched vortices are anti-correlated and 
their total vorticity vanishes. 

We recognize from Eq. (^), that only neutral con- 
figurations of vortices will have finite energy. Hence the 
partition sum can be restricted to neutral configurations 
(i-*^- X]r ■'^R. — 0) like in the pure case. Every such neu- 
tral configuration can be imagined as a superposition of 
dipoles of a vortex {N — 1) and an antivortex {N = — 1). 
In analogy to electrostatics, vortices of the same sign re- 
pel each other and vortices of different sign attract each 
other. 

Since the interaction is logarithmic, vortices build a 
Coulomb gas. The interaction is identical among the 
thermal and quenched component. The Fourier represen- 
tation ( pa| ) of the vortex Hamiltonian makes evident, that 
thermal vortices try to compensate quenched vortices as 
well as possible. A perfect compensation is prevented by 
the fact, that A'^r is restricted to integer values, whereas 
Qr is a continuous random variable. 

In Eq. ( |3b|) we introduced the vortex core energy 
E — T:"fK, where 7 takes the value 7 w 1.6 on a square 
lattice. It emerges from the lattice Greens function [cf. 
Eq. (4.13a) of Ref. |l^. Although 7 has a unique value 
for the lattice model, it will occasionally be considered 
as an independent parameter which allows to control the 
density of vortices irrespectively of K. In a dilute vor- 
tex system with A^r = 0,±1 only, E is comparable to a 
chemical potential. 

In a region of parameters, where vortices are negligi- 
ble, fluctuations of the original angle variables are given 
only by spin wave&-||-|They lead to a decay of the spin 
correlation functionllEl 



r(r) = (cos[0(r) - ^(O)]) - r"'' (5) 

with exponent rj — {1/K + a)/2TT. If the correlation de- 
cays algebraically even in the presence of vortices, the 
system has quasi long-range order. 

As vortices are excited thermally, they give rise to ad- 
ditional fluctuations in the angles, the correlation will 
decay faster. This effect can be expressed by a renormal- 
ized rj, which becomes scale dependent. If vortices lead to 
a divergence of the renormalized rj on large scales, quasi 
long-range order is destroyed. However, vortex fluctu- 
ations are not only responsible for such a quantitative 
effect, they also drive the phase transition into the high- 
temperature phas%, where the spin-spin correlation de- 
cays exponentially.Q 
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Let us estimate, to what extent disorder can favor the 
excitation of vortices. In the pure case Q = 0, the ground 
state of (^) is vortex-free. The density of vortices van- 
ishes for vanishing temperature, since the creation of 
quantized vortices always costs finite energy. Therefore 
the pure system has a phase with quasi long-range order 
at low temperatures. 

A simple argument shows, that the disordered sys- 
tem has_ft5£fin at zero temperature a finite density of 
vorticesOtj we consider two positions at distance R 
and determine the probability of finding there a vortex- 
antivortex pair. The energy of such a dipole is given by 
Udipoic{R) ~ 27rJlniZ, where we neglect the core energy 
for large R. Disorder gives and additional energy contri- 
bution Vdipoic(^). At T = 0, the dipole will be present, 
if its total energy J7dipoio(^) + Vdipoic(-R) is negative. The 
original Gaussian distribution of the random phase-shifts 
results in a Gaussian distribution of Vdipoic(J?) with vari- 
ance A2(i?) := V^~JRj ~ 47rcrj2 Ini?. Thus the prob- 
ability 



It is convenient to introduce the ordinary, disconnected, 
and connected correlation, 



P{R) 



dV 



27r2 In R 



x/2TA(i?) 



(6) 



of finding a pair at T = is finite. In other words, this 
shows that in the presence of disorder the ground state 
is no longer vortex-free and that renormalization effects 
can be important even at arbitrarily small temperatures. 



III. SCREENING 

The goal of the present work is therefore a calculation 
of such renormalization effects in the presence of ran- 
dom phase-shifts. We use the conceptual most simple 
approach, a self-consistent linear response approach-Jor 
the pure Coulomb gas, this approach is equivalenfJiZi to 
the real-space renormalization group of Kosterlitz.B We 
have checked that this equivalence still holds in the pres- 
ence of disorder. 

We now define vortex renormalization effects by 
macroscopic properties of our system. As a probe we 
introduce additional test vortices into the system. They 
experience a screened interaction and the screened back- 
ground potential acting on them. Let us denote the 
unscreened interaction between vortices t/k = Att^K-^. 
Due to quenched disorder vortices Q, thermal vortices 
N are subject to a Gaussian background potential Vk = 

The screened 



^^kQk with variance V-k^ 
background potential and interaction are identified from 
the contributions to the free energy of test vortices, which 
are of first and second order in their vorticitylla (see also 
Appendix ^ for some intermediate steps) 

c/k(gk- (A^k)), 



= C/k-C/k^^co„(k). 



(7a) 
(7b) 



C(R) ~ (N-EcNo), 
Cdi«(R) (ArKAoI, 
C'con(R') := C'(R) — Cdii 



(R) 



(8a) 
(8b) 
(8c) 



For a specific realization of disorder, the screened inter- 
action in principle depends explicitly on the coordinates 
of both interacting vortices, it is not translation invari- 
ant. Although we should in principle consider screening 
in the particular disorder environment, we approximately 
evaluate only the average screening effect by taking the 
disorder average, which restores the translation symme- 
try of the interaction. 

From the large-scale behavior of the screened interac- 
tion and the variance of the screened potential one can 
then identify screened parameters K^'^'^ and a^^'^ as (for 
the precise definition see Appendix ^) 

/CxD 
dRR^Ccon{R), (9a) 

/GO 
dRR^Cdis{R) + 

/oc 
dRR^Ccon{R)- (9b) 

Although we do not yet know the correlation functions, 
we may establish from these expressions criteria for the 
existence of an ordered phase. In this phase, we expect 
vortices to give rise only to finite screening effects, i.e. 
the screened parameters must have finite values. This 
requires lim^j^oo -R^C'con(-R) < oo and a similar condi- 
tion for the disconnected correlation. 

In the limit of large core energy, one can calculate the 
correlations to leading order by considering only a sin- 
gle vortex dipole. In th absence of disorder, one finds 
Cco„(R) = C(R) ~ and Cdis(R) = 0. The con- 

dition for the ordered phase then simply reads K > 2/7r, 
in agreement with the condition that the free energy a 
single vortex has to be positive.Q A similar argument can 
be constructed for the disordered system at T = 0, where 
one can estimate C(i?) ~ p{R) ~ using Eq. (||). 

Then the condition for order reads a < 7r/8. Jitis ar- 
gument already disproves previous predictions,EJl3'Ll that 
infinitesimal cr > would destroy order at T = 0. 

For finite core energy, screening effects have to be taken 
into account quantitatively for the calculation of the cor- 
relations. For this purpose we develoi)|-a self-consistent 
scheme in analogy to the pure case.D'ta We introduce 
scale-dependent variables K = K{1) and a = (t(Z), which 
include screening by dipoles of radius 1 < i? < e' only. 
ThenEq. (|) can be cast in differential form: 



^i^-i = -27rVCeo„(i? = e'), 
|a = -27rVCdis(i? = e'). 



(10a) 
(10b) 
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Both equations combine to a simple flow equation of the 
spin-spin correlation exponent, -^rj = — 7r^e^'C(i? — e'). 

Eqs. ( [lO| ) suggests a simple "two-component picture" 
illustrating the screening effects of vortices. This picture 
virtually separates dipoles of thermal vortices N into a 
frozen and a polarizable component. The frozen compo- 
nent gives rise to the Cdis- It is not polarizable and does 
not contribute to screening of the interaction. However, 
since it is frozen, it behaves like the quenched disorder 
background and leads to an effective disorder strength 
represented by the screened a. The polarizable compo- 
nent is the source of Ccon and leads to the screening of 
K. Both components contribute equally to a renormal- 
ization of the spin exponent. This virtual separation into 
two components must not be taken literally. It is impos- 
sible to assign a particular dipolc uniquely to one of the 
two components. The two-component picture is analo- 
gous to the two-fluid model of superfluidity, which also 
must not be taken literally in the sense that a particular 
atom would be either superfluid or normalfluid. 

In order to take advantage from the differential screen- 
ing equations, we will calculate the correlations at i? = e' 
in a self-consistent way which includes screening effects 
by smaller dipoles. At each differential step we will also 
perform a simultaneous rescaling R e^'^'R. Then we 
can interpret the differential equations as renormaliza- 
tion group flow equations. 

For the further investigation-of the model, we apply 
the standard replica techniqueEj to the partition func- 
tion (^. We are going to rewrite screening in the repli- 
cated system in order to expose the consistency between 
screening in the unreplicated and replicated system. 

The Coulomb gas ( |3b| ) is replicated n times and after 
disorder averaging one obtains the effective Hamiltonian 

R ab 

R^R' a,b 

We introduce the couphngs K"'' = KS"'' - K, K ^ 
aK^/{l + naK) and the core energy E"''' = tt^/K'''' = 
ES"^ ~ E, where E = ttjK and E = iT'jk. On the 
partition sum the neutrality condition 

ENr = (12) 

R 

is imposed for every replica. 

In the replicated system, screening effects can be calcu- 
lated in the same linear response scheme as before. The 
leading term of the free energy of test replica- vortices is 
of second order in the infinitesimal test vorticity. From 
this order, we derive a differential screening 

j^K"'' = 27r3e4' ^ K'"'C'"^{R = e})K'^'' (13) 

cd 



in analogy to Eq. (|lC|) . We introduced the replica corre- 
lation C^R) := {N^N^)n: where (. . .)„ denotes a ther- 
mal average in the replicated system. Since the initial 
interactions are replica-symmetric, i.e. the coupling has 
the form K'^'' = K5°'^ — K, we use the replica symmet- 
ric ansatz C"'' = don'?"'' + C^^-^ for correlations. The 
relation to the correlations in the unreplicated system is 
provided by 

Cx(R)= lim4")(R) (14) 

for all three types X (ordinary, disconnected, and con- 
nected) of correlations, if we identify C^"^(R) := C""(R) 
consistently. 

Exploiting the replica symmetry of correlations, Eq. 
( |T3| ) decays into flow equations for K and of K, which 
read in terms of K and tr = K /{K"^ — nKK): 

j^K-^ = 47r3y2^„, (15a) 

j^a - Air'yl,. (15b) 

In order to approach the usual notation of renormal- 
ization groups, we define effective fugacities 

yj,^yj,{l):=-^e*'cP{R=l) (16) 

again for all three types of rescaled correlations. They are 
related by = j/^j^ -I- y^^^. In these definitions we took 
already account of rescaling of lengths, which generates 
an additional flow of the core energies, 

j^E = ttK, (17a) 

4;E = nK. (17b) 
dl 

Eqs. ( p^ together with ( p7| ) constitute the main renor- 
malization group flow equations. From relations (^4|) we 
recognize, that the flow equations of K and a of the un- 
replicated system coincide with those of the replicated 
system in the limit n ^ 0. In the next section these 
equations will be completed by a flow equation of the 
free energy, which however does not feed into the other 
flow equations. 



IV. SELF-CONSISTENT CLOSURE 

The flow equations (|l^) together with ([l7| ) can not yet 
be evaluated since they are not closed. In this section 
we are going to achieve this closure by expressing the 
effective fugacities in terms of K, a, E and E. 

In the replica language, the effect of disorder is encoded 
in additional interactions between vortices and we first 
have to discuss their nature. An inspection of Eq. ( p4| ) 
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shows, that vortices of opposite sign in the same repHca 
interact with a potential 2Tr{K — K)lnR. For n > 1 
they always attract each other. These dipoles are stable 
against dissociation due to thermal, fluctuations only for 
(K—K) > 2/7rMEl In previous workElfl the stability of pre- 
cisely such dipoles was used as criterion to determine the 
boundary of the ordered phase. However, for n = and 
low temperature or strong disorder, K — K = K ~ aK^ 
becomes negative (repulsion!) and this criterion becomes 
questionable. The increasing instability af these dipoles 
gives rise to a reentrant phase boundary.EI 

However, it is not sufficient to consider interactions 
only within replicas. Independently on temperature and 
replica number n, vortices of the same sign but in differ- 
ent replicas attract each other with a potential 2'kK In R. 
In the limit n ~ Q the coupling between different replicas 
reads K = aK^. For low temperature or strong disor- 
der, the interaction between different replicas becomes 
stronger to the same extent as the intra-replica cou- 
pling becomes weaker and correlations have to be calcu- 
lated taking into account the competition between intra- 
and intcr-rcplica interactions. The importance of intefc 
replica interactions has been recognized by Korshunov.Q 

The essence of the replica problem is a suitable treat- 
ment of the inter-replica interaction for integer n > 1 
and to construct an analytic continuation which allows 
to take then the limit n ^ in a proper way. In the fol- 
lowing this is done in a way different from Korshunov's 
way, leading to opposite conclusions. On the technical 
level, this step is the essential progress of the present 
work. 



A. Physics of n > 1 

As necessary prerequisite for the analytic continuation 
n — > 0, we must discuss vortex fluctuations in the repli- 
cated system for integer n > 1. This happens in some 
detail since our final conclusions contradict Refs. 

We first check, that the ground state of the repli- 
cated Hamiltonian is the state without vortices. For this 
purpose we examine the Hamiltonian in Fourier space, 

Hn /kEah(27rV'^^)^"''^-k^k- There we easily rec- 
ognize, that every creation of vortices costs energy, be- 
cause K°'^ is positive definite: K°'^ has one eigenvalue 
Ki = K - k = K{1 + [n- l)aK)/{l + naK) > and 
n— 1 eigenvalues K2,...,n = K > 0. Since the ground state 
has no vortices, we are allowed to use for low tempera- 
tures a low density expansion. This is in contrast to the 
unreplicated system, where the ground state has a finite 
density of vortices, as discussed above. 

For the following it is convenient to introduce the no- 
tion of a "replica- vortex" (compare Fig. Q): in a state 
with vorticity {iV^, . . . , N-^ at position R, we say that 
a "replica-vortex" of type v = {N^, . . . , N^} is at po- 
sition R. (Upper indices at N are not exponents but 
replica indices!) Since the core energy contains a term 



^^Rai^R.)"^^ which suppresses large |iV"|, we may re- 
strict our consideration for large or 7 to states with 
— 0, ±1: one vortex-free state and 3" — 1 different 
types of replica- vortices, antivortices included. 

Because of the neutrality condition (|l^) the elemen- 
tary excitations in the replicated system are again vor- 
tex dipoles. In the dilute limit (large 7), some type 
of such replica-vortices will form with their antivortices 
those dipoles, which are most instable against thermal 
fluctuation and thus destroy order. To find out this 
type, we have to discuss the energy of such dipoles. The 
core energy of replica- vortex 1/ = {N^, . . . , iV"} de- 
pends only on the total vortex number mg :— X]a(^°)^ 
and the number mi :— ^ (mo -I- iV°) of vortices with 
N"- — +1. The same holds for the coupling K,y, which 
gives the strength of the logarithmic interaction between 
a replica-vortex of type i' and its antivortex: 

= nioE - (mo - 2mifE, (18a) 
= m^K ~ (mo - 2mi fk. (18b) 

Due to replica symmetry, there is a certain degeneracy 
between different types which we include when we 
speak of the class {mQ, mi) of replica vortices. 

The existence of order requires stability of all types of 
dipoles against thermal dissociation. For the dilute sys- 
tem this is guaranteed by miui^^o > 27r according ta 
the simple stability criterion of Kosterlitz and Thouless.cl 
The least stable class with the weakest coupling can be 
determined explicitly: it is given by replica-vortices of 
class (mo — 1, mi = 0, 1) with K^, = K — K in a, high 
temperature region (n-\-l)K < K , and by replica-vortices 
of classes (mo = n, toi — 0, n) with Ki, — nK — n?K in 
a low temperature region {n + 1)K > K. The resulting 
phase diagram for n — 1,2, 4, 8, 16 is shown in Fig. |^. 

The fact, that the class of the least stable replica- 
vortices for 71 > 1 depends on the physical parameters, 
can be interpreted as an indication for a phase transition 
even for n — 0. We take it as a warning, that we must 
not focus on a single class of vortices_The origin for the 
failure of the early replica approachealrO is just the focus 
on classes (mo = l,mi = 0,1) only. We even do not 
restrict ourselves to the two classes which are dominant 
for n > 1, since it turns out a posteriori, that also other 
classes contribute to the replica limit n ^ . 



B. Towards the rephca limit n ^ 

Now we calculate the correlations and the effective fu- 
gacities taking into account contributions by all types of 
replica- vortices. For a low density of replica-dipoles (re- 
alized for large 7 or at low temperature) we can neglect 
the interaction between different replica-dipoles ("inde- 
pendent dipole approximation" ) . This means that we can 
calculate approximate correlations in the partition sum 
containing only up to one dipole of replica- vortices. 
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R„ — Ti 



(19) 



The fugacity of replica- vortex 1/ is defined hy yi, := 
g2;-_B„ rp]-^g factor 1/2 is present, since v runs over vor- 
tices and their antivortices and we should count every 
reahzation of a dipole only once. In Eq. (|l9|) integration 
is restricted to |R+ — R^;] > 1 for alH > 0, since lengths 
are rescaled. 

Before we turn to the calculation of correlations, we 
wish to make sure that calculations in this ensemble make 
sense in the limit n — > 0. For this purpose we have to 
convince ourselves, that Z„ ^ 1. This is equivalent to 
the condition, that the free energy per replica has a fi- 
nite limit. Then also the renormalization group flow of 
the free energy per replica, which follows from the con- 
tribution of dipoles of radius 1 < i? < e"^' to Z„, is finite. 
The flow reads 



4^ = 2^-.lVy^ 



dl 



(20) 



where J- denotes the free energy per unit volume, per 
replica, and divided by temperature. The first term on 
the right-hand side originates from rescaling. 

From Eq. ( pO| ) one can derive (see Appendix that 
a restriction to a single class of replica-vortices leads to 
an inconsistency with the replica trick, namely the diver- 
gence of for n ^ 0. Therefore we retain all classes of 
vortices. 

Before we turn to the correlations, we evaluate the sim- 
pler sum, which is the source for the fiow of free energy: 



e^' {i^- 1}, 



(21) 



where we introduced a "dummy" Gaussian random vari- 
able A with ^ = and ^ \/2 and a "shell partition 
sum" z = z{A, E,E) := 1 + + z_ with weights 



z± 



-2(E±A-\/2E) 



(22) 



Physically, the "dummy" Gaussian random variable 
A represents nothing but the disorder component which 
couples to the dipoles with radius in the shell being in- 
tegrated out. Since disorder is uncorrelated on different 
length scales, we can average over disorder on a given 
scale right when we consider screening by dipoles with a 
radius of just the same scale. For this reason we can undo 
the replica trick in the flow equations in the very same 
way as we introduced it initially in the full unrenormal- 
ized system. In Eq. (^ the replica number n became 



an explicit variable and we may now perform the analytic 
continuation n — > 0, 



dl 



„4/ 



In ; 



(23) 



As shown in Appendix this yields a well defined free 
energy per replica only because we have retained all 
classes of replica- vortices. This means, that types of fluc- 
tuations, which seem to be irrelevant for n > 1 (in our 
case: energetically expensive replica-dipoles) can become 
important in the replica limit. 

In a similar way we can proceed to determine the effec- 
tive fugacities. In the independent dipole approximation 
( |l9| ) the rescaled correlation functions read (for i? > 1) 



(24) 



We ignored normalization by the factor since this 

factor becomes unity for n = 0. With the help of a gen- 
erating function 



„4i 



e 



C({^}) 

we find (a 7^ 6) 



-2E ^^(jV»)^+^^(2^V2^+»?'' 



(25) 



drfdrf^ 



drfdrf 



jj=0 



jj=0 



„4i 



(^++^-) 
„2-n 



(26a) 
(26b) 



The i?-dependence of correlations can be restored by the 
substitutions E ^ E + ttK In R and E ^ E + ttK In R. 
As before, n became an explicit variable such that we can 
send n ^ 0. The effective fugacities then read 







z+ 


+ Z 




1 + z 


+ + 


1 


ylis = 


¥' 


(if 


f ^ 






z+ 




Vcon 




Z+ + 




+ 4z+z_ 


(1 + 


z+ 


+ Z-)2 



(27a) 

(27b) 
(27c) 



Now we have achieved our goal of expressing the ef- 
fective fugacities of the disordered system (n ~ 0) as 
functions of E and E. In principle these fugacities could 
also depend on K and u. Such a dependence is absent in 
the independent dipole approximation used above: cor- 
relations are determined neglecting interactions between 
dipoles, and the energy of single dipoles of radius unity 
does not depend on K and a. 

Equations (15), (p^), and ( |27| ) form a closed set of 
equations. The flow of free energy (|23) does not feed 
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back into the other flow equations. Appeahng to the two- 
component picture, we may call "density of dipoles" , 
/ := e^'lnz/2y^ "reduced free energy per dipole", p :— 
Vconlv^ "fraction of polarizable dipoles" , and q :— 1—p — 
Vdis/y'^ "fraction of frozen dipoles" . These quantities are 
given by Eq. ( p7| ) as functions of E and E. For the con- 
venience of the reader we summarize the flow equations 
with the two-component parameters: 





—T 




dl 


d 














d 








dl" 




'^E 








dl^ 




d - 








dl^ 



2T-2^fy\ 


(28a) 


iTT^py^, 


(28b) 


An^y^, 


(28c) 


ttK, 


(28d) 




(28e) 



As initial values for K and a one has to use the unrenor- 
malized values, which also enter E = ttjK, E — tt^cjK'^ 
and thereby y, f, p, and q. Initially JF(Z = 0) = since 
no fluctuations are included. 



C. Asymptotic approximation 

In the present form the flow equations are not very 
convenient since the dependence of the effective fugaci- 
ties and the two-component parameters on E and E is 
still quite intricate. Therefore we perform an additional 
approximation which is valid on large length scales. 

First the integrals over A (which are to be performed 
as disorder averages in Eqs. (pTf)) are split into inter- 
vals where one of the contributions to the shell partition 
sum dominates. Then the normalizing denominators in 
Eqs. ( ^ ) can be expanded with respect to the smaller 
contributions. The resulting infinite series can be inte- 
grated over A analytically term by term. From these 
series one then can extract easily the leading terms for 
large / after approximating E « ttKI and E « iraK^l. 
Since -^E = nK and -^E = iraK^, these approximations 
are asymptotically correct provided K and a converge to 
finite values. 

For convenience we introduce an effective temperature 
variable 



1 



2aK 



(29) 



The asymptotic approximation yields for the flow of the 
reduced free energy 



,2.F-yif5^e-('^/2-)'+4' (r<l) 
For the fugacities we find analogously 



dT 
Id 



(30) 



-2ir/S'(l-l/2T)/+4i 

'(■n /2<y)l+il 



V 277^ sill 



(r>l) 
(r < 1) ' 



Ucon 



■ sin(7rr) 

, l) (r > 2) 

(-<2) 

g-27rK(l-l/2T)i+4i ^ 2) 

V 2n^l sin(7rT) ^ ^ 



(31a) 
(31b) 
(31c) 



Although we were starting from unique expressions, 
the asymptotic behavior differs for high temperature or 
weak disorder and for low temperature or strong disor- 
der. There are three regimes of parameters, namely r > 2 
(called regime lA), 2 > r > 1 (called regime IB), and 
T < 1 (called regime II). Their physical properties will 
be analyzed in the next section. 

Approaching the separatrices r = 1 or r = 2 using the 
expressions for the high- and low-temperature regime, 
the expressions should coincide. This is true for their ex- 
ponents. The comparison of the prefactors is not mean- 
ingful, since in the low-temperature expression the di- 
verging factor l/sin(7rr) is multiplied with the asymp- 
totically vanishing factor 1 / Vl. 

Remarkably, for zero temperature (r — 0) we find 
y'^ = y^on = R^P, which means that the probability P, 
Eq. (^, coincides with the ordinary and connected cor- 
relation. 

For the analysis of the flow equations, we wish to use 
the quantities of the two-component picture, as intro- 
duced in Eqs. (^). The reduced free energy per dipole, 
fraction of polarizable dipoles, and fraction of frozen 
dipoles converge asymptotically to 



p := 



q := 



41 In 2 




2y^ 




Ucon ^ 








Vdis _ f 


' 


y^ ^ \ 





1 (r > 1) 

1/t (r < 1) 

(r>l) 
(r<l) ' 

(r>l) 

T (t < 1) ■ 



(32a) 
(32b) 
(32c) 



Here we have neglected contributions which vanish on 
large scales. For example q is not strictly zero but van- 
ishes exponentially with I — > oo for ct > and r > 1. 
Remarkably, the asymptotic two-component parameters 
depend on K and a only through the effective tempera- 
ture T. 

In order to eliminate E and E completely from the 
asymptotic flow equations, we determine the flow equa- 
tion for the fugacity: 



dl^ 



2y + ttK^ + ttK^ 



(4 



dE dE 
2ttK{1 - aK))y^ 



y 



(r>l) 
(r<l) 



(33) 



To close this subsection of the paper, we wish to point 
out, that the flow equations (28) for K and (33) for y 
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with expression (321:) for the polarizability coincide with 
those of Nattermann et aE3 apart from factors smaller 
than 2. This small deviation is due to a rougher asymp- 
totic approximation in Ref. 



V. RESULTS 

In the previous sections the technical part of deriving 
the flow equations has been achieved. The most impor- 
tant flow equations are Eqs. (p8|a-c) . Originally, the two- 
component parameters /, p, and q have been functions 
of E and E. In the asymptotic approximation E and 
E have been eliminated, the two-component parameters 
are given by Eq. ( |3^ ) as simple functions of r = l/2aK, 
and the fugacity flows according to Eq. (p3|). These flow 
equations can be summarized by 



d 2 


= (4 - 2'KT*K{\ - aT*K))y^ 


(34a) 




= 4:n^T*y^ 


(34b) 


d 

dl" 


= in^il ~ T*)y'^ 


(34c) 


d ^ 
—T 

dl 

r 


= 2T - 27r— y2 
min(l, r) 


(34d) 
(34e) 



Thus we are in position to discuss the physical implica- 
tions of these flow equations. 



A. Qualitative aspects 

To start with, we address qualitative aspects and show 
that the flow equations satisfy some fundamental physi- 
cal criteria. 

According to Eqs. (|l^) and (^^ all correlations Cx are 
negative, since all fugacities are positive y\ > 0. This is 
consistent with a screening that weakens the attractive 
interaction between vortices of the same sign and weak- 
ens the repulsive interaction between vortices of the same 
sign. 

As a consequence, the two-component parameters are 
also non-negative. Our asymptotic expressions show, 
that the fraction of polarizable dipoles p = t* vanishes 
linearly at small temperature, whereas the fraction of 
frozen dipoles q = 1 — t* becomes unity. 

Screening reduces the coupling J — KT since ^-f^ < 0, 
and it increases the strength of disorder, > 0. In com- 
parison to the pure case the screening of J is reduced by 
the factor < p < 1, the fraction of polarizable dipoles. 
Whereas previous theoretical approaches did not yield a 
renormalization of cr, Eq. (|3^c) manifests that our renor- 
malization scheme leads to an increase of a. The free en- 
ergy density is reduced by screening, see Eq. (||d). Due 



to disorder the contribution of dipoles is modified by a 
factor / = 1/r* > 1. 

The increase of a seems to contradict the tendency of 
thermal vortices N to neutralize quenched vortices Q. 
This neutralization is expected from Eq. How- 
ever, the screened a is not defined by the fluctuations 
of Q— (N). Rather, it was introduced via the fluctuation 
of the screened disorder potential V. The fluctuations 
of V are proportional to aK^. The expected neutraliza- 
tion shows up as -^crK^ < 0. This holds even though 
4tT > since the reduction of K dominates the increase 

dl — 

of cr, which follows from 2aKp > q that is easily proven 
with the asymptotic expressions. 

The most important quantity to analyze is the flow 
of entropy. From spin-glass models one knows, that a 
replica-symmetric theory might lead to a negative en- 
tropy. This unphysical feature then indicates that the 
symmetric approach is insufficient and one has to allow 
for a breaking of replica symmetry. We now show that 
the present replica-symmetric theory does not suffer from 
this problem. Therefore we believe that a breaking of 
replica symmetry is not necessary. 

When the entropy S = —-^TT is derived from the 
reduced free energy, the unrenormalized values of J and 
a have to be kept fixed. According to the flow equa- 
tion for T ^ a non-negative entropy is guaranteed by 
ln{Tfy'^) > 0. The demonstration of this inequality is 
complicated by the temperature dependence, which is im- 
plicit to the screened J and a for ^ > 0. If one ignores this 
implicit dependence, one easily verifies hi(Tfy^) > 0, 
since ^Tf > and ^y"^ > hold for r > 1 and r < 1 
separately. The implicit temperature dependencies can 
be taken into account by showing ^^ln(r/j/^) > 0. 
This is verified in the leading order, where one has 
^ HTfy^) = (4 - 27rif (I - aK)) + 0{y^) for r > 1. If 
one would use the same expression for t < 1, one would 
obtain a negative entropy flow. However, for r < I we 
find jj HTfy^) = (4 - ^) 0(y2). Therefore the en- 
tropy flow is non-negative even at very low temperatures. 



B. Ordered phase 

A numerical integration of the flow equations (even 
without asymptotic approximations) shows, that below 
some transition line Lc in the (T/ J, i7)-plane the renor- 
malization flow indeed does converge for Z — > cx) to fixed 
points with Koo, Coo, and yx = 0. This is true not only in 
the limit of large 7 (small fugacity), but also for arbitrar- 
ily small 7. Only the extent of the ordered phase shrinks 
with increasing 7. This implies in particular, that the 
original XY- model described by 7 « 1.6 has an ordered 
phase. Above Lc all quantities diverge under renormal- 
ization and the present renormalization group leaves its 
range of validity, which is limited to small fugacities al- 
ready due to the independent dipole approximation. 

The properties of the ordered phase are evident in 



8 



terms of the renormahzed parameters Koo and ctqc • Due 
to screening they are related to the bare parameters by 
Koo < K and (Too > cr, where the inequalities become 
equalities for vanishing fugacity within the ordered phase. 

In the ordered phase, the irrelevance of vortices re- 
quires that the fugacities vanish on large scales. There- 
fore the flow equation (^Ja) of fugacity must have a nega- 
tive eigenvalue. This is true at temperatures and disorder 
strength below L^, which is located at 




(2/7r <Koo< 4/^) 
(i^oc > 4/^) 



(35) 



in terms of the renormalized parameters. It is depicted 
in Fig. H as function of the renormalized parameters and 
of the unrenormalized parameters for ^ — 1.6. 

The most striking feature of the phase boundary is the 
absence of reentrapce. Reentrance was predicted by ear- 
lier approaches. H'LI Absence of reentrance without a loss 
of the orderedjphase was found only recently in improved 
approaches .EJI13 

As a function of the renormalized parameters, the 
phase boundary crcio(^oo) = 7r/8 has a horizontal part 
at large i^oo, i-e. at low temperatures. In terms of 
bare parameters, the line can be parametrized by a 
function a^K) or Kdcr). Due to the renormalization 
of disorder strength, ^ct > and due to the reduction 
of the coupling -^K < the end points are found at 
ac{T = 0) < 7r/8 and Kc{a = 0) > 2/7r and crc(T) has 
no longer a horizontal part for finite fugacity. 

Concerning the asymptotic behavior of the vortex 
correlation functions, one can distinguish three regions 
within the ordered phase. They are separated by a line 
Lci, fixed by 2croo^oo = 1 (equivalent to Too = 1), and 
a line Lc2, fixed by 4CToo-fi'oo = 1 (equivalent to Too = 2), 
see Fig. |[ In the (J/T, cr)-diagram these lines emerge 
from the origin. They are straight lines in the limit of 
zero fugacity and bent upwards for finite fugacity. 

The correlations are retrieved from the fugacities, Eq. 
( ^ ) , by omitting rescaling and inserting I — \nR. 
Thereby we find (i? > 1) 



C(R) 



-2R 



-2 



27r^ In R sin(7rTo, 



Ccon(R-) ~ \ ^ I Oo 



27r^lni? sin(7rroo) 



(36a) 
(36b) 



for Ttx) ^ 1 and 



Cdis (R-) 



-2 



TTToojl - Too) 



27r2 1ni? sin(7rroo) 



for Too ^ 2. The limiting values for the fraction of polar- 
izable and frozen dipoles are related to the correlations 

by 



lim 

R^oo 

lim 

R — *oo 



Ccon (R) 
Cdis (R) 

C(R) ■ 



(37a) 
(37b) 



Naturally the question arises, whether these regimes 
are separated by a true thermodynamic transition at Ld 
and Lc2- As we have observed, that asymptotic form 
of the flow equations and of the correlations changes its 
dependence on a and K non-analytically at t = 1 and 
T = 2. However, a phase transition in the strict sense 
requires a non-analytic dependence of the integrated free 
energy on the unrenormalized parameters a and K . We 
were not able to prove such a singularity at Lci and Lc2, 
which we therefore consider as crossover lines only. How- 
ever, a true transition is expected at Lc just as in the 
pure case. 

Region II can be distinguished qualitatively from the 
other regions. One can consider qoo as order parameter 
for the low-temperature regime II. It measures the rela- 
tive strength of the disconnected vortex density correla- 
tion to the usual correlation on large scales. In the def- 
inition of this order parameter the normalization of the 
disconnected correlation by the usual correlation is nec- 
essary, since true long-range order is absent. The glassy 
order reflected by the Edwards- Anderson like correlation 
can therefore be only quasi long-ranged, too. Neverthe- 
less one may compare qoo to the Edwards- Anderson order 
parameter in spin-glasses and we interpret it as order pa- 
rameter for the glassiness of in the vortex system. 

The regions lA and IB differ only by the analytic de- 
pendence of the connected correlation on Koa and CToo. 
The correlation under consideration do not provide an 
order parameter comparabe to q for the line Lc2- 



C. Universality 

The universality class of the pure XY-model is defined 
by the behavior of the critical exponents 77 = 1/4 and 
(5 = 15 at the transition and the exponential divergence of 
the correlation length when the transition is approached 
from above Jj 

In the presence of disorder, the spin correlation decays 
with exponenlE 



Vo 



0/2^. 



(38) 



(36c) 



Along the liiie Lc the exponent takes the value 7700 = 1/4 
for a = 0,0 which decreases with tempetature until it 
reaches the value rjao = 1/16 for T = 0. til In the dis- 
ordered phase, where correlations decay exponentially, 
Tjoa = 00 formally. Thus the value of rjc^ at the transi- 
tion becomes non-universal, i.e. it is disorder dependent. 
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Only for cr = or T = the jump of 7700 is universal 
(from 1/4 or 1/16 to 00 respectively) in the sense that it 
does not depend on the bare coupling J or bare fugacity 
y. For cr > or T > the jump lies inbetween and will 
depend on microscopic details. 

In the present work we have not included a magnetic 
field. Therefore we can not determine 5 in the presence 
of disorder. 

We now attempt to analyze how the correlation length 
^+ diverges when the transition line is approached from 
above. Even though the flow equations become invalid 
for large fugacity, KosterlitzU has argued that ^+ is given 
by the length scale where the fugacity becomes of or- 
der unity. In the pure case, where the flow equations 
reduce to jjK'^ = Air^y^ and = (4 _ 2nK)y'^ he 

found i+ - exp(6/\/i^-i - K^^) with a non-universal 
constant b. In the zero-temperature limit the flow equa- 
tions read ^cr = An'^y'^ and Jjj/^ = (4 — j^)y'^- By a 
replacement a — > l/AK they can be mapped onto the 
flow equations of the pure case. Therefore we conclude, 
that the zero temperature transition has 



^+ ~ exp(fe/\/CT - cTc) 



(39) 



with some other non-universal b. For finite temperature 
and finite disorder, our flow equations are too compli- 
cated for being integrated up analytically. However, we 
naturally expect that the functional form of the diver- 
gence is preserved and that only the argument of the 
square-root has to be replaced by the distance to the 
transition line, say min(3.^_^^)gi^[(A'^^ — ^cT^)^ + ^ 



VI. DISCUSSION AND CONCLUSIONS 

The present theory is based on a selfconsistcnt linear 
screening approach, which is solved in a differential form. 
This differential form has been presented in the renormal- 
ization group language. In fact, the very same differential 
equations can be directly derived by a real space renoij=. 
malization group, generalizing the scheme of KosterlitzEl 
to the replicated system. In this scheme the flow equa- 
tions obtain contributions from integrating out all classes 
of replica dipoles with the smallest (cut-off) diameter. 

The controversy among earlier theoretical treatments 
of the disordered XY-model requires to put the present 
work into perspective with these earlier approaches. 

When Cardy and Ostlund (CO) considered the more 
complicated problem of the XY-model in the presence 
of random flelds,El they used a replica-symmetric ap- 
proach to derive renormalization group flow equations, 
which generated also random bond disorder. Rubinstein, 
Shraiman, and Nelsoru (RSN) focused on the random 
bond problem, which they analyzed without using the 
replica formalism. However, the resulting flow equations 
were a subset of those of CO. Their starting point was 



essentially Eq. (10a), which they evaluated only to the 
leading order in the fugacity of vortices in the unrepli- 
cated system. Consequently, they neglected the discon- 
nected correlation and could not find a fiow of cr. Their 
fiow equations for K and coincide with our asymptotic 
fiow equations for r > 1. However, they used them also 
for r < 1, which gave rise to reentrance and a negative 
fiow of entropy. 

The present approach is in principle equivalent to CO 
and RSN. The essential difference lies in the calculation of 
the correlations, or in the contributions which are taken 
into account when the dipoles of smallest distance are 
integrated out. CO and RSN thereby made too crude 
approximations, which treat the disorder effectively as 
annealed. This yields correct results only for high tem- 
peratures r > 1 but fails for t < 1. 

Technically, they considered only contributions of 
replica-vortices of class (mo — 1,toi — 0, 1). We pointed 
out, that already for replica numbers n > 1 this restric- 
tion is insufficient at low temperatures. Therefore we 
included all classes of replica- vortices. The additional 
vortex classes turned out to be important for r < 1 and 
remove reentrance and the negative flow of entropy. 

From the disappearence of reentrance after the inclu- 
sion of additional vortex classes we have to conclude, 
that these additional fluctuations in the replicated sys- 
tem (71 > 1) reduce the effect of thermal fluctuations in 
the original unreplicated system {n = 0). This strange 
feature reminds of the fact, that in variational replica ap- 
proaches one has to maximize and not to minimize thfij 
free energy for n = Oin order to flnd the physical state.c2l 

When Korshunovu considered the problem in the 
replica representation, he attempted to identify the tran- 
sition in the limit of zero fugacity with the dissociation of 
the least stable class of dipoles. The A^-complexes in his 
notation are identical to dipoles formed by replica vor- 
tices of class (mo, mi = 0, mo), where one has to identify 
= mo. They indeed drive the transition for n > 1. 
From my point of view, there are problems in his way of 
taking the limit n — > 0, which cause the conclusion, that 
the ordered phase should not exist at all. One problem 
is, that it does not make sense to consider replica vor- 
tices which occupies mo > 1 replicas even for a number 
of replicas n < 1, since (n — mo) (mo — 1) > should be 
preserved during the analytic continuation. 

In addition, the consideration of such restricted classes 
of dipoles is not sufficient for an analytic continuation 
n ^ 0. Even then, if one takes into account all classes of 
replica- vortices, it is the most stable and not the least sta- 
ble class, which drives the transition in the limit n — > 0. 
This is again a strange feature of replica theory. The 
arguments for these statements are detailed in Appendix 

In Appendix |y we further show, that Lc can be de- 
rived already from considering a single vortex in a finite 
system. After replication its free energy can be calcu- 
lated by taking into account certain types of fluctuations, 
which are again parametrized by some m in the range 
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(n — m) (to — 1) > 0. The transition follows from 
a maximization of the free energy with respect to this 
parameter. 

For several disordered systems the occurence of reen- 
trance and negative contributions to the flow of entropy 
have been shown to be artefacts of a replica-symmetric 
theoretical approach, -which could be removed by break- 
ing replica symmetry.tj For the present model these arte- 
facts have emerged in the early replica-symmetric ap- 
proaches, indeed. Here we have shown that these prob- 
lems can be cured by including the relevant fluctuations 
within a symmetric theory. Hence a breaking of replica 
symmetry is not necessary. 

The recent apntnaches of Nattermann, Scheidl, Kor- 
shunov, and LiliilrJ and also of Cha and Fertigi3 found 
the existence of the ordered phase and the absence of 
reentrance without using the replica formalism. Their 
results agree with the present theory at least in the limit 
of zero fugacity, for finite fugacity there are quantitative 
deviations. In a very recent real-space renormalization 
group approach in the unreplicated system, Tang£J has 
derived flow equations identical to Eqs. (Q) and achieves 
even quantitative agreement with the present results. 

A qualitative progress of the present theory is the 
inclusion of a renormalization of a. This renormaliza- 
tion is driven by the disconnected vortex density corre- 
lation, which in earlier approaches was neglected from 
the very beginning or found to be zero due to later ap- 
proximations. Although a is increased under renormal- 
ization, this effect does not destroy the ordered phase. 
The flow of a is essential for the universality type of 
the zero-temperature transition, which deviates from the 
Kosterlitz-Thouless transition in the pure system only by 
the value of the critical exponents but not in the expo- 
nential divergence of the correlation length. 

Due to the renormlaization of a one can exclude a hor- 
izontal part of Lc at low temperatures even for finite 
fugacity. Such a horizontal part was shown to exist by 
Ozeki and NishimoriEl for a wide class of disordered sys- 
tems with a certain gauge symmetry. But the XY-model 
under consideration does not possess this symmetry. 

To sum up, is has been shown that the two-dimensional 
XY-model with random phase-shifts on bonds has a 
phase with quasi long-range order. It has exhibits a zero 
temperature transition at a critical strength of disorder. 
Near the transition the correlation length diverges expo- 
nentially as in the pure case, but the jump of the spin 
correlation exponent is of different magnitude. 

Within the ordered phase there is a glassy low- 
temperature region. In the absence of true long-range or- 
der, we identified the ratio of the Edwards- Anderson cor- 
relation to the usual vortex-density correlation on large 
length scales as order parameter for glassiness. 

We used a replica-symmetric renormalization group 
approach. The reentrance and negative entropy flow of 
earlier approaches of that kind could be cured by tak- 
ing into account additional fluctuations. No need for a 
breaking of replica symmetry was found. 



Several generalizations of the present work are of inter- 
est. In a forthcoming publicationp3 the present theory 
will be extended to include (random) magnetic fields on 
the ordered phase of the XY-model. The properties of 
the system in the disordered phase are beyond the range 
of our perturbative renormalization approach. In partic- 
ular the question about possible glassy properties in this 
phase deserves further research. In the limit cr — > oo, the 
present model becomes the so-called |gauge glass model, 
which also exhibits a glass transition.E3 
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APPENDIX A: LINEAR RESPONSE 



Since Eqs. (10) are the basis of the present analysis, it 
seems worthwhile to present some intermediate steps of 
their derivation. We start by including a test vorticity q 
into the effective vortex Hamiltonian (pa|), 

= i^C/k|iVk + 9k-Qk|'. (Al) 
The free energy of the system with test vortices reads 

^[{g}] = -ln^e-«^. (A2) 

{N} 

The screened potential and interaction are defined by 

5T[{q}] 



5q-i 



9=0 



[/^"■(5(k + k') 



S-^T[{q}] 



5q--kSq-i,' 



(A3a) 
(A3b) 



q=0 



and result in Eqs. (^. The screened parameters are 
identified via 



:= hm — C/^"-, 



(A4a) 
(A4b) 



Fourier transformation then leads to Eqs. (H). On the 
way one needs the probability distribution of quenched 
vortices, 



P[{Q}] oc exp 



1 



(A5) 
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and one has to integrate partially 



gk(iV_k) - J V{Q}P[{Q}]Qk{N^k) 



47r2 

U2 



-ai^C7con(k). 



(A6) 



The resulting screening formulae are also valid in the 
case, where the interaction and disorder correlation on 
small scales have additional contributions due to E and 
i5, which flow under renormalization independently of K 
and a. 

It is important to realize, that a^'^'^ could alternatively 

be defined by ^{Q-k - {N^k)){Qk - (iVk)) for small 
fc, which also might be plausible at first sight but gives 
different results as soon as K^" ^ K . Then, however, it 
would be necessary to distinguish the screened interac- 
tions among thermal vortices N from the screened inter- 
actions between thermal vortices N and quenched vor- 
tices Q. This complication is avoided by the definition 
actually used, which is based on the correlation of V^'^^ . 
Only in this way screening effects can be captured by a 
renormalization of the parameters K and cr, which were 
present already in the unrenormalized model. 



APPENDIX B: EXTREMAL PRINCIPLE 

In this appendix we address the question, whether the 
correct vortex flow equations can be found, if one con- 
siders only the contributions of a single class of replica 
vortices to the partition sum and the correlations. 

For n > 1 one naturally expects, that fluctuations of 
the vortices with the lowest energy cost give the most 
important contribution to renormalization. Since the 
interaction between different replicas favors vortices of 
the same sign, optimal replica vortices will be of class 
{mQ,mi = 0) or (mo,mi = toq); which we call simply 
m = mo vortices. According to Eq. ([l8|), their core 
energy is 



Em = niE - 



m'^E. 



(Bl) 



Since E and E are positive, the dipoles with minimum 
energy are always given by m = 1 for {n '\- l)E < E and 
by m = n for {n + 1)E > E. They determine the phase 
diagram of Fig. I The most stable pairs are given by 
TO* = E/2E (or, for n > 1, by one of the two nearest 
integers to m*). 

Now we consider screening effects by m-vortices only. 
We would like to point out, that the restriction to this 
class of vortices does not break replica symmetry, since 
we include all 2(^) possibilities to place such vortices in 
n replicas. Then the flow equations read 



al n\m I 



11^^ ^" tfconJ 



(B2a) 
(B2b) 
(B2c) 



with fugacity y„ 



and effective fugacities 



(™_\)2/m, Vlis = il-DvL and y^^^ = {mJi)y?n- They 
can be related by the identities (J^) = (".l') = 

and ("-?,) = i;2-l("-i). Thereby we find 

n — 1 Vm— 1/ ' \m — 2/ n— 1 Vm— 1/ 

the two-component parameters and flow of y 



TO 



n-1 ' 

TO — 1 



(B3a) 
(B3b) 
(B3c) 
(B3d) 



where Km '■= rnK — rn^K. 

To simplify the discussion of the replica-limit n = 0, 
we use the asymptotic approximation E « -kKI and 
E K, TTdK^l. One should consider to in the range 
< TO < 1, which follows from an analytic continuation 
of the condition (n — to) (to — 1) > 0. 

One easily convinces oneself in the case without disor- 
der [E — Q), that one obtains unphysical results, if one 
looks for TO which minimizes Km'- this would be to = 
and, the systen would be unstable to the creation and 
dissociation of dipoles {Eq = Ki^ — Q). However, if we 
look for the maximum^ which lies in the pure case always 
at TO = 1, we find a positive core energy Ei and the 
correct phase boundary at Ki = K = 2/Tr. 

Therefore let us try to look for the most stable type 
of dipoles also in the presence of disorder. It is given by 
TO* = T := \/2aK. Recalling the constraint < to* < 1 
in the limit of n = 0, we have to distinguish a high- and 
low-temperature regime with 

m* = 1; Km* ^ K{1-1/2t) for t > 1, (B4a) 
m* = t; Km' = l/4cr for r < 1. (B4b) 

If the core energy is large and renormalization is negligi- 
ble, the criterion Km* > 2/7r for the dissociation of the 
most stable dipoles then yields the correct transition line 
Lc as in Eq. (^5|) . Surprisingly, even the two-component 
parameters / — 1/r , p = and g = 1 — r for r < 1 ac- 
cording to Eq. (B3) are in perfect agreement with (^). 

So, was it essentially unnecessary to consider all classes 
of vortices? No! Even though the fiow equations are 
correct, we must know the initial values to integrate 
them up, in particular we must know the unrenormal- 
ized value of y. The problem is, that for n — > we have 
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2/^ = (™._\)2/m* ""It ^ ym": which diverges - n ^ in 
the low temperature regime with < m* < 1. This pre- 
tends erroneously the relevance of vortices. In particular, 
there is an inco nsis tency with the replica trick, since ac- 
cording to Eq. ( |B2| a) the free energy per replica diverges 

These prob- 



like 



cx n 



1 / n— 1 \ ^ sin(7rm*) 

m* Vm*— 1/ m*7rn 

lems are removed by including the contributions of all 
classes of vortices. 



T = {ttK{1~ <tK) -2}\nL, 



(C7) 



whereas for iraK^ > 2 the minimum is taken at 
^Jlj-KoK"^ < 1 and we obtain 



J- — ttK 



Ini. 



(C8) 



APPENDIX C: SINGLE VORTEX ARGUMENT 

It is instructive to show, that the phase boundary Lc 
in the limit of zero fugacity can be derived by consider- 
ing a single vostex. Similar to the reasoning of Kosterlitz 
and Thoulessp we examine the possibility of finding sin- 
gle free vortices in a system of finite linear dimension L 
by determining its free energy. 

The energy of a vortex is composed of its self-energy 
irJlnL and of a potential energy V due to disorder. The 
disorder-averaged free energy reads 



T = TTK\nL + 5T, 



ST = -In / d^re-P^i'') 



lim — In Z" 

n^Q n 



Z" = / dVi . . . dV„ X 



(CI) 
(C2) 

(C3) 



We use y^(ra) w InL and find the leading power 

of Z" in L by considering all possible partitionings V of 
replicas into groupsa 

(C4) 

Here \V\ denotes the number of groups, and (5^^ = 1 or 
if replicas a and b are in the same or in different groups. 
If we consider only the possibility, where n/m groups of 
m replicas are formed, we find 



(C5) 



The restriction to these partitionings looks at first sight 
like a one-step replica symmetry breaking, but it is not. 
If one partitioning is included, also all other partition- 
ings, wich are related by a permutation of replicas, are 
also included. 

In the analytic continuation n —^ the group param- 
eter has to be taken in the range < m < 1 and the 
maximization is replaced by a minimization. Then 

ST = - min \2/m + TraK^m] In L. (C6) 

0<m<l '- ^ 

Finally, for iraK'^ < 2 the minimum is taken at m — 1 
and we obtain 



The phase boundary is there, where the free energy 
of the single vortex changes sign. Eq. ( |C7|) gives the 
high-temperature part of Lc, whereas Eq. ( jCq ) gives the 
low-temperature part of Lc in agreement with Eq. (|3^). 
Note, that the boundary iraK^ = 2 separating the va- 
lidity of Eqs. (C7) and ( |C8| ) does not coincide with Ld 
of Lc2- This, however, cannot be expected since within 
the ordered phase fluctuations are dominated by vortex 
dipoles and not by single vortices. The appearence of the 
separatrix naK^ = 2 is specific for the physics of single 
vortices. 
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FIG. 2. Phase diagram for n = 1,2,4,8,16 in the limit 
7 — » oo. Below the transition line all types of dipoles are 
stable against thermal dissociation, above the line there is 
at least one instable type. At high temperatures (above the 
shoulder of the transition line) , the transition is driven by the 
dissociation of replica dipoles (mo = l,mi = 0, 1), whereas it 
is driven by dipoles (mo = n, mi = 0, n) at low temperatures 
(below the shoulder of the transition line). 
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FIG. 1. Schematic picture of a dipole of replica-vortices. 
Every horizontal line represents a two-dimensional system, 
which is replicated n — 8 times in this example. An up/down 
arrow represents a vorticity A*'r = ±1. The dipole con- 
sists of a vortex of type f — {0, 1, 1, 0, —1, 1, 0, —1} and class 
(mo = 5, mi = 3) and its antivortex at positions R = R*. 
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FIG. 3. Phase diagram for 7 = 1.6 (thick lines) and 7 = 10 
(thin lines) . The transition Lc (continuous lines) separates the 
low-tompcraturo phase with quasi long-range order from the 
disordered high temperature phase. The crossover line Lci 
(long dashed lines) separates the glassy region II from the 
non-glassy region I. On the crossover line Lc2 (short dashed 
lines) separates regimes lA and IB within phase I. There the 
disconnected vortex density correlation function changes its 
asymptotic behavior. 
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